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Observations indicate galaxies are distributed in a filament-dominated web-like 
structure. Numerical experiments at high and low redshift of viable structure 
formation theories also show filament-dominance. We present a simple quanti- 
tative explanation of why this is so, showing that the final-state web is actually 
present in embryonic form in the overdensity pattern of the initial fiuctuations, 
with nonlinear dynamics just sharpening the image. The web is largely defined 
by the position and primordial tidal fields of rare events in the medium, with 
the strongest filaments between nearby clusters whose tidal tensors are nearly 
aligned. Applications of the cosmic web theory to observations include probing 
cluster-cluster bridges by weak gravitational lensing. X-rays, and the Sunyaev- 
Zeldovich eff"ect and probing high redshift galaxy-galaxy bridges by low column 
density Lyman alpha absorption lines. 

The large scale structure we observe in the distribution of galaxies is thought to be a 
consequence of the gravitational instability of a medium that was relatively smooth in the early 
Universe, with the nonlinear patterns we see having evolved from small density fluctuations 
superposed on a homogeneous and isotropic background. The basic elements in the pattern 
are clusters, rich and poor, filaments and sheets, connected in a network, with voids largely 
devoid of galaxies encompassing the volume in between.^ A^-body simulations evolving the 
best-motivated classes of initial fiuctuations into the nonlinear regime also show a network of 
filaments, a result that was visually evident for many years, ^'"^ but not really understood. A 
sample "cosmic web" A^-body pattern is shown in Figure I, for group and cluster formation in a 
"cold dark matter" model Universe: as the density threshold drops from high values the regions 
which first emerge are clusters, then arms stretching from the clusters, which ultimately join 
to form the predominantly filamentary network: the first pattern to percolate is filamentary 
and it is that which the eye picks out. The basic filamentary pattern is seen on smaller scales 
at higher redshift in simulations, with bright galaxies playing the role of clusters. There was a 
suspicion that either the structure was an artifact of simulations or would need fully nonlinear 
dynamics to understand.'* 

Two classical stories have largely defined how people have thought about structure in the 
Universe, the (Russian) pancaking picture^'^ and the (Western) hierarchical clustering picture. 
The Zeldovich theory^ applied to models with fiuctuation power on cluster scales but none 
on galactic scales. It led to a highly infiuential description of the medium, based upon the 
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catastrophe theory of the caustics in the density field®: planar pancakes form first, these 
drain into filaments which in turn drain into clusters, with the entirety forming a cellular 
network of sheets. The rival hierarchical clustering theory applied to models with initial 
density fiuctuation power decreasing with increasing scale, with a characteristic mass scale 
for nonlinear objects increasing with time.^''' The main paradigm used here was a merging- 
halos picture, where halos are associated with peaks in the 'background' linear (Gaussian- 
distributed) density field, ^° i.e., of the rare events in the medium. 

The nonlinear dynamical evolution of an initially cold medium is described by a map, 
x(r,t) = r — s(r,t), from Lagrangian (initial state) space r to Eulerian (final state) space x. 
In the early linear phase, the displacement field s is small and the map is one-to-one, but 
eventually when particle orbits cross in the medium the map becomes multivalued. We show 
here that the filamentary web is a consequence of the distribution and spatial coherence of 
the strain field in the medium. To establish this story, we shall split the displacement field 
s = Sf, -\- Sf into a background field S;, smoothed over a large scale Rj, and a fiuctuating field 
Sj built from shorter wavelength components. We also characterize the filtering by the linear 
rms density fiuctuation amplitude smoothed on that scale crp[Rb^t). The conceptual picture is 
one of a smooth stately large-scale single-stream fiow S;, upon which is superposed a relatively 
chaotic nonlinear short-distance fiow Sj in which orbit crossings abound; we must restrict 
(Tp(i?6,t) to be ^ 1 for this picture to be valid. The linear approximation (in Lagrangian 
space), Sb[r^t) = D[t)sb[r)^ where D[t) describes the linear growing mode of fiuctuations,'* 
defines the "truncated" Zeldovich map, which by itself gives a reasonable description of the 
overall pattern of A^-body simulations. 

The 'background' strain field is defined by eb^ij[r) = — (VjSfcj -|- VjS6i)/2. In linear theory, 
the shear tensor is just its time derivative, fifc^ij, is proportional to the gradient of the 
smoothed gravitational potential and hence the smoothed peculiar tidal field acting on a patch 
of matter is proportional to e^^ij. The strain tensor has eigenvalues —Xyj^ which we order by 
Xv3 ^ Xv2 ^ ^vi- We also express^^ the Xyj in terms of the height of the density fiuctuation Pj, = 
[Xyi-\-Xy2-\-Xy3) / (Tp^ the shcar ellipticity e^, = (A^,3 — A^,i)/(2z/6(Tp) and prolaticity = (A^,i-|-A^,3 — 
2At,2)/(2z^6crp); thus e^, > and —Cy < py < ey. In this notation, extreme sheet-like structures 
would have py ^ Cy, extreme filaments would have py ~ —Cy. Both are statistically unlikely: 
in the conditional distribution^^ P({A^,j } Iz/^), {cy Iz/^) 0.54;/^" is significantly nonzero while 
{Pv\'^^b) is zero (and the dispersion is ~ 0.2u^ for both); and we find^^ regions with py ^ 
and large Cy favour filaments over sheets. The conditional distribution P{{XyEj}\Szc) of the 
Zeldovich-mapped shear eigenvalues,^"^ XyEj = D[t)Xyj[[l — D[t)Xyj]~^ ^ given the Zeldovich- 
mapped overdensity at the percolation threshold, Szc (where (1 -|- 6z) = 11(1 -|- A^,^;^)), is 
qualitatively similar to P[{Xyj}\pb)^ and further favours filaments. 

Fig. 1(b) shows the structure of the isodensity surface of the Zeldovich background field 
map at the percolation threshold. Even with this approximate dynamics, the web defined by 
filaments that was evident in Eig. 1(a) is clearly seen: sheets are not the dominant elements. 
In the truncated map with full nonlinear rather than Zeldovich dynamics, the cluster regions 
occupy less Eulerian volume which leads to an even more enhanced filamentary character than 
in Eig. 1(b). And the prominent filaments in the final state are there in fattened form in the 
initial state, as Eigure 3 shows. 

To quantitatively probe the coherent structures in the initial conditions, we turn to a 
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new language: conditional multiple-point correlation functions in Lagrangian space — i.e., 
statistically-averaged density and displacement fields subject to various constraints on the 
shear at multiple points. Consider first the behaviour of an ambient patch about a single 
random point in the medium at which we specify the shear parameters Cy and py to be the 
typical ones corresponding to a given height Pj, in the Lagrangian space: we find the typical 
density contours are more filamentary than sheet-like even before Zeldovich mapping and the 
density coherence scale along the f-direction, which is the filamentary axis, is greatly enhanced. 

The remarkable size of the filaments is not derivable from the one-point constraint. For 
this, we turn to correlations constrained by at least two rare peak-patches in the medium. 
Examples are shown in Figure 2. The peak-patch theory^^ shows that it is statistically likely 
that, given a specific orientation of the shear tensor for a peak-patch, neighbouring peak- 
patches will be preferentially along this f-axis and have shear tensors aligned with it^"*. We 
allow the separation to vary as well as the shear orientation. What we see is that the strong 
correlation bridge between two clusters gradually weakens as the separation increases, or as 
the shear tensor orientations become misaligned. The reason for this phenomenon is that the 
high degree of constructive interference of the density waves required to make the rare peak- 
patches, and to preferentially orient them along the 1-axis, leads to a slower decoherence along 
the 1-axis than along the others, and thus a higher density. 

So we paint the following picture how clusters and their filamentary bridges weave a web: 
Clusters nearer than the mean cluster separation can have strong correlation bridges. There are 
many such pairs since clusters are statistically biased. ^° The bridge breaks if the separation is 
too large, isolating the clusters from each other — unless there is a cluster in between to which 
both have extended their filamentary bridges. To make this into a network, consider laying 
down the rare events in the medium according to the clustering pattern of peak-patches which 
become clusters when they evolve dynamically. The correlation bridges arch from cluster 
to cluster in much of the domain, and these dynamically evolve to filaments, creating the 
network and containing the bulk of the mass. The typical scale of the segment (bridge) of the 
filamentary network is ~ 30h~^ Mpc. The order in which the physically significant structures 
arise is basically the inverse of that in the classical pancake picture: first high density peaks, 
then filaments between them, and possibly afterwards the walls, defined as the rest of the mass 
between voids. 

An important byproduct of this picture is our ability to sequentially capture more and 
more large-scale structural features by imposing constraints from a rank ordered list of peak- 
patches. It gives us a method for reconstructing initial states using the rarest peaks (clusters) 
and, through a background-field map, the main features of the final states. The excellent 
reconstruction shown in Figure 3 develops as peaks are sequentially added, as is shown in 
Figure 4. The anisotropic peak-strain constraint is essential for success in reconstruction. 

The list of peak-patches (and void-patches) is therefore a powerful way to rank order and 
maximally compress the information stored in the initial conditions, showing what is essential 
to define structures on the basis of a modest set of local measurements. The characterization 
is spatially localized: only the rare events nearby define the region. And a modest extrap- 
olation of the work reported here confirms the conjecture of ref. [11]: within a cluster-scale 
peak-patch, the rarest group-scale peak-patches defined at moderate redshift will define the 
dominant evolutionary characteristics (the merger history), and one has to work one's way 
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down the group-patch hst only if one wants finer and finer dynamical detail. This has im- 
portant applications to the construction of constrained realizations for nonlinear numerical 
simulations. 

The web pattern depends upon the shape of the initial state power spectrum, 
which is usually parameterized by a local power law index f^p,e//(^) = din {\6L[k)\'^) / din k. 
In our figures, we used the standard CDM spectrum, which has rip^eff ~ —1.2 on cluster 
scales. When we steepen the spectrum to rip^eff ^ 0, the clusters become less clustered and 
the coherence of the web is lost, and although some filaments remain they are weaker and 
shorter. When we fiatten it to rip^eff ^ ~2, the clusters are more clustered, the coherence is 
more pronounced and the filaments are both strengthened and widened. We also see stronger 
membranes in the regions between the filaments when a number of clusters are close together. 
These are sheet-like structures, but are not classical pancakes. 

There are several immediate applications of our theory to observations. We suggest super- 
clusters are predominantly cluster-cluster bridges; no image from the cosmology of the 80's 
was as powerful as the CfA picture^^ of the Coma cluster and its environs, with a great arm 
of galaxies spanning Coma and Abell 1367. "Great Wall" features would be interpreted as the 
membrane-like webbing in between the filaments, and would be more pronounced the fiatter 
the spectrum; e.g., theories which match cluster observations have rip^eff ~ —1.6 on cluster 
scales, hence on average would be more membrane-like than CDM would be. 

We predict that the strongest filaments exist between highly clustered and aligned clusters. 
This should be especially notable in the mass distribution around systems of clusters recon- 
structed with weak gravitational lensing,^'' which uses the projected strain (shear) field. It 
also may affect the strategy of the search for X-ray gas^® or the Sunyaev-Zeldovich effect in 
filaments. Most theories predict density fiuctuation spectra are fiatter at galactic than cluster 
scales, with rip^eff ^ ~2 appropriate for redshifts z ^ 3. The neutral hydrogen absorption by 
the filaments existing then would be observed as part of the quasar Lyman alpha forest, with 
column densities Nhi ^ lO^'^'^cm"^ and a characteristic (comoving) filamentary size roughly 
the distance between dwarf galaxies, ~ lh~^Mpc. As well, we expect filaments of dwarf 
galaxies between normal galaxies at similar redshifts. 
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Figure Captions 

Figure 1. (a) A^-body simulation (courtesy of A. Klypin) of a cold dark matter (CDM) 
model with 128"^ particles in a comoving periodic boundary volume of (50 li~^Mpc)"^ evolved 
to the rms density fluctuation amplitude erg = 0.7 (where erg = dp at a top-hat filtering scale 
Rb = 8h~^Mpc, where h parameterizes the Hubble constant, Hq = fOOhkm s~^ Mpc~ ). f in 
20 particles are plotted. The pattern of connected filaments joining clusters evident here is 
quite general for viable hierarchical theories that explain the galaxy clustering data. 

(b) Truncated Zeldovich map for the (Lagrangian space) initial conditions that led to the 
Eulerian space pattern of (a). The initial state was smoothed on a Gaussian scale of Rj, = 
5h~^Mpc, corresponding to a linear amplitude of (Tp = 0.54. The density contour level shown, 
Szc = 1, corresponds to the percolation threshold. At high threshold [6z ^ 4), the contours 
are isolated around rare massive clusters, at medium threshold [6z ~ 2), arms radiate from the 
clusters, and there are some web-like membranes close to the clusters; by Szc = 1, the filaments 
connect the clusters. Making Rj, smaller shows more detailed structure, and the filaments 
can break up into beads upon them (the eye still catches the coherence), but the picture is 
otherwise similar. As we see here, the Zeldovich map does not imply that pancakes associated 
with caustics are the dominant overdensity enhancement structures in the medium for (Tp{Rb) = 
(9(f/2). (We have also checked whether the contours of XvE?, bear much relationship to the 
observed structures, as might be expected in the classical pancake scenario: they do not.) 
It is the rare positive and negative peaks that condition the overall pattern formation in the 
medium. We identify the filaments that are so visually evident as "correlation bridges" between 
the collapsed peak-patches, forming a web, but no classical pancakes. 

Figure 2. Zeldovich-mapped correlation functions constrained by shear at two peaks show the 
effect of changing separation on the correlation bridge and changing shear-tensor orientation. 
What we see is that the strong correlation bridge between two clusters gradually weakens as 
the separation increases, breaking by 50h~^Mpc, or as the shear tensor orientations become 
misaligned. For the left panels, shear tensors are aligned, and the separations are 40h~^Mpc 
(top), 30h~^Mpc (middle) and 50h~^Mpc (bottom, with the broken bridge). For the right 
panels, the separation is 40h~^Mpc and the misalignment angles relative to the line joining 
the peaks are: ±30° (top), ±20° (middle) and ±90° (bottom), when the bridge is fully broken. 
(At this 8zc = 1 level with ±45° orientation the bridge is also broken, but joins at a lower 
threshold.) Mean-shear parameters appropriate for Abell-cluster peak-patches were adopted: 
V = = 0.18 and = 0. (The average separation of clusters of this richness is 50 h~^Mpc 
and their correlation length is ~ 20h~^Mpc; i.e., clusters come in patches. ^^) To avoid too 
much interior dynamics in the patches, the truncated Zeldovich map was used with Gaussian 
smoothing scale i?6 = 5h~^Mpc {(Tp = 0.54), the same as for the peaks (whose top hat scale is 
Rpj^ 10h~^Mpc). Note that in Lagrangian space the length of the filament bridging the two 
peak regions is similar to the diameter of the peak patches. The contour shown is at 8zc = 1, 
about where percolation occurs in the A^-body computation at this (Tp smoothing. 
Figure 3. The upper left panel shows the linear CDM density field 6l{'^) in Lagrangian space 
for the (50h~^Mpc)"^ box of Fig. 1. It has been smoothed on a Gaussian scale i?6 = 3.5h~^Mpc, 
corresponding to (Tp = 0.65. The threshold chosen is 8^ = l(Tp because the linear density field 
percolates at this one-sigma level. The upper right panel, a repeat of Fig. 1(b), is the Zeldovich- 
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map of the smoothed initial conditions into Eulerian space, for a contour threshold bz = 2, 
just above where percolation occurs. The final state filaments are present in the initial state. 
The lower left panel shows the linear density correlation function ((5/,(r) |20 pks) constrained 
by the top 20 peak-patches identified in the initial conditions, and the lower right shows its 
Zeldovich map. 

This mean field reconstruction of bi^ given "observations" of the values of the shear tensor 
and the displacement at N peak-patch points Fpfc averaged over the peak-patch size 
((5/,|A^pks). This is a straightforward application of a well-known theorem for a constrained 
Gaussian field, ref. [7], Appendix D. This gives 9A^ constraints. The requirement that the 
gradient of the density field vanishes at Fpfc (the extremum constraint) adds a further 3A^ 
constraints, but is unimportant in practice. Although the positions are not counted explicitly 
as constraints, changing them radically alters the way the mean fields look. The greater N is, 
the smaller are the fiuctuations about the mean. 

The patches are indicated by blackened ellipsoids of overall size proportional to R,p}^ and 
shape defined by the shear tensor orientation, with the shortest axis corresponding to the 
highest shear eigenvalue, the longest to the least. [To find peak-patches^^, we first identify 
peaks in the smoothed linear density field above a threshold bi^^ = 1-686 (the linear overdensity 
needed for a spherically symmetric "top hat" to collapse to a point); the smoothing is done 
on a hierarchy of Gaussian filter scales, here, 25 filters ranging from 1 to 4h~^Mpc. This of 
course greatly overcounts sites of local collapses — a very strong version of the "clouds-in- 
cloud" problem that is present even for single-filter peaks, as noted in ref. [10]. About each of 
these points, we calculate the radius R,p}^ at which the interior region is predicted to have just 
collapsed in all three dimensions at the redshift in question. At a higher redshift, the radius 
would be smaller, and at a lower redshift larger. For this we use the homogeneous ellipsoid 
modeP^ for interior dynamics which utilizes both the interior tidal field and a linear-evolution 
approximation to the exterior tidal field acting on the patch. This requires measuring the 
strain (or shear) tensor within this radius (top-hat filtered). The technique^^ also gives the 
local displacement field averaged over the Lagrangian volume within. In the next stage of 
the peak-patch procedure, smaller scale peaks whose centres lie within larger scale peaks are 
removed from the list. The net list is rank-ordered in mass. We find about 400 peak-patches 
in our (50h~^Mpc)"^ volume, with the top 20 ranging in size R,p}^ from 4.5 to 9.2h~^Mpc, with 
masses from lO^'*-^ to lO^^'^M®.] 

Figure 4. The top three panels show mean fields with an increasing number of peak-patch 
and void-patch constraints imposed. Left shows Lagrangian space, right shows the result after 
Zeldovich-mapping into Eulerian space. The top panels use the top 5 peaks, the second use 
the top 5 voids (treated as "inverse-peaks") as well as the top 5 peaks, and the third use the 
top 10 peaks. Adding 5 voids does not do as well as adding 5 more peaks. The fourth set, 
using constraints on the peak height and the velocity (or displacement field) only, shows how 
important the anisotropic shear constraints are for defining the structure. The peaks used are 
shown as solid ellipsoids, the voids as open ellipsoids, scaled in size oc with shape defined 
by the shear tensor, shortest axis corresponding to the highest shear eigenvalue, longest to the 
least. At the thresholds 6l = 1^ 6z = 2 of Eig. 3 some strong bridges seen in the 20 peak 
reconstruction, ((5/,|20 pks), are not as evident in the ((5/,|10 pks)-field, but emerge at lower 
thresholds; 6l = 0.5 and 6z = 1.3 are shown here. This figure shows that it is better to 
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use the top 2N peaks rather than N voids and N peaks if we are focusing only on regions 
of higher density contrast — which is observationally most relevant for most astrophysical 
observables. The situation would be reversed if it were voids we were concentrating on. It 
is rather interesting to see how the void-patch constraints create high mean field regions in 
between them, just where less rare peak-patches reside. It is just not as precise a way to get 
that structure for the same computational effort. 
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